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^ l' We explore simulations on periodic lattices in the Tomboulis SO(3) X Z{2) for- 

5^ ' mulation. We measure gauge invariant vortex counters for "thin", "thick" and 

^^ "hybrid" vortex sheets in order to tag Wilson loops by the occurance of gauge in- 

, variant vortices linking them. We also measure projection vortex counters defined 

" ' ■ in the maximal center gauge for comparison. 

►^ . 1 Introduction 

^^ ' Lattice QCD continues to maintain an important role in the search for the 

f— ^ . physics of color confinement. The lattice regulator maintains gauge invariance 

^\ ' at all costs. Dynamical variables are group elements rather than elements of 

^^ . a Lie algebra. As a consequence many of the topological features that are 

^P ' prominent candidates for elucidating the physics of confinement have natural 

^ ! lattice definitions. These include U(l), Z(N), SU(N)/Z(N) monopole loops, 

"t^ ' Dirac sheets, Z(N) and SU(N) vortex sheets etc. These objects are often abun- 

'— H , dant in U(l) and SU(N) lattice gauge theories. They become singular only as 

O ' one approaches the continuum limit. 

D . Yaffeld, TomboulisEl and Kovacs and TomboulisQ have developed a formu- 

"^ ' lation of SU(N) gauge theory that is manifestly SU(N)/Z(N) invariant. In this 

^ . formulation, center elements, Z(N), multiplying each link leave the action and 

k> ' measure invariant. New Z(N) variables, defined on plaquettes, cr(p), carry the 

'V^ . Z(N) degrees of freedom. This formalism, equivalent to the standard SU(N) 

C^ ' form, allows an elegant topological classification of the SU(N)/Z(N) and Z(N) 

vortex configurations occuring on the lattice. See also Refs. 4-7. 

This talk is a report on our paper 13, in which we address the issues of 
simulations in the Tomboulis variables on a periodic lattice. Many results 
have followed from this formulation without doing simulations El. As a first 
calculation, we tag Wilson loop measurements by the occurance of vortices 
linking the loop using a number of different vortex counters. We restrict our 
attention to SU(2). 

We also measure the P (projection) vortex counter in the original SU(2) 
formulation following Refs. 10 - 18 for comparison. Projection vortices arise 
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in a Z(2) gauge theory derived from the original SU(2) theory by going to the 
maximal center gauge and then replacing links by sgn{tr{Ux^p,)). The projected 
theory has "thin" Z(2) vortices defined on one lattice spacing. They have been 
found to be correlated with center vortices and therefore a measurement of P 
vortices is a predictor of them. 

2 Z z (2) xsu {2) /z {2) on atoms 

The Wilson form of the partition function, 

Z = J [dUib)]exp f/3^ itr[C/(ap)] j , (1) 

can be cast into a form in which the variables transform under Z{2) and 
SU{2)/Z{2) = S0{3) separatelyH. 

Z = J [dUm J2 CWipHp)] exp hj2l\trmdp)]\a{p)\ , (2) 



ct(p) 



where 



C[cT{p)fjip)] =J2l[ Xrip) HpMp)) n ^ (^(^^) ■ (3) 

r(p) P b 

The derivation makes use of the invariance of the Haar measure under Z(2) 
transformations U{b) — > U{h)^{h) of the links. Summing over 7(6), the link 
matrices become SU{2) representatives of the invariance group SU{2)/Z{2) = 
SO{i). The Z{2) dependence is transferred to new variables defined on pla- 
quettes o'(p). The variable r]{p) is the sign of the plaquette tr[C/(9p)] = 
|tr[t/(i9p)]| X 'q{p). 

C[cr(p)?7(p)] is a constraint, equal to either a positive constant or zero. 
The Z{2) character, Xa{P) = ilj where a — ±1 labels the representation 
and /3 = ±1 is the argument. x+(+) = x+(-) = X-(+) = -X-(-) = +1- 
The delta function constrains the r summation to contribute only if an even 
number of negative t plaquettes are in the coboundary of each link. 

For free boundary conditions □ 

C[q;(p)] = TT<5 (a(9c)) X constant 

c 

This requires that there be an even number of a{p) — a{p) x 77 (p) plaquettes 
on the faces of all elementary cubes. 



3 Zero weight configurations on the torus 

This constraint C[a(p)] is further restrictive for periodic boundary conditions. 
It gives zero for a vortex configuration of negative plaquettes on a cocfosed 
surface wrapping around the torus. This means that one can not excite a 
configuration corresponding to the topologicaUy stable anti-periodic boundary 
conditions and vice versa. _ 

To see this we use a resuh from Ref. U Note that the {r} configurations form 
a group with muhiphcation rule {tiT2){p) — ti{p)t2{p). Using the invariance 
of the group sum over r we can write 

chp)] = e n ^r(p) («(p)) n ^ (^(^^)) ' 

t{p) P b 

t(p) P b 

= {to, a) X C[a{p)l 

where (r, a) = -^ Yin Xt{p) ('^(p)) r • If we can find a group element tq for which 
(to, a) ^ 1 thenC[a] =0. 

Take the example of a coclosed stack of a = —1 (1,2) plaquettes at lo- 
cations (0,0,^,77) and a closed tiled surface of r = —1 (1,2) plaquettes at 
locations (/i, i^, 0, 0). They have one plaquette in common, giving one factor of 
— 1 in (to, a), and hence C[a] — 0. 

As a second example let us consider how C[a] gives the cube constraint. 
Let us suppose that a particular cube violates this constraint with an odd 
number of faces with a = —1. Then take a configuration tq which takes 
values —1 on all 6 faces of this particular cube. This is a closed surface and 
therefore satisfies the constraints imposed on r. Again IIpXtoCp) (q^(-P)) = ^1 
and therefore C[a] = 0. 

4 a{p) configuration space 

We are interested in simulating in the variables {U{b),a{p)}. Allowed con- 
figurations a{p) are defined indirectly by the constraint, function C[a\. The 
corresponding simulation would be cumbersome to implement. 

We proposedEl a constructive definition of allowable a{p) configurations by 
building them up from "star transformations" , i.e. correlated sign flips of the 
a plaquettes occurring in the co-boundary of each link. We showedls that the 
constrained updates of six plaquettes reaches all allowed {a} configurations. 



and in fact is identical to the above definition. In this section we summarize 
this result. 

Before discussing the a configurations let us first describe the link updates. 
This is a straightforward generalization of the link updates for SU{2). The 
proposed change in a link might change the sign of the 77 plaquettes in the 
co-boundary of the link. If one of these changes sign, we need to flip the sign 
of the corresponding a plaquette so that the a configuration is unchanged. 
Then the Monte Carlo step is essentially the same as for the SU(2) update. 

Next consider the above mentioned star transformations. Our proposed 
update is to flip the sign of the six a plaquettes forming the co-boundary of 
the links. 

It is easy to see that both these update steps will preserve the local con- 
straints, i.e. the cube constraints. 

Next examine the-global constraints. Consider the operator constructed 
out of /i, v plaquettes H 

■^^''^ = n ^(P)^(P) - v{S^.jy)cr{S^,^^) = ±1, 

where 5*^,,^ is a whole tiled /i, v plane. Afp,^i, = ± for an even/odd number of 
vortices of stacked fi, v plaquettes wrapping around the orthogonal ^,77 direc- 
tions of the torus. 

We start with A/'p,,y = +1 in all 6 planes. It is easy to see that our update 
algorithm preserves A/"^,!/. 

Let us define the relevant sets of configurations more carefully. 






reC p 



5^(r,a) =. J2(a,r). (4) 



TEC TEC 



The bracket, {ol,t) = ±1, is negative if and only if there are an odd number 
of plaquettes for which both a{p) = — 1 and t[p) = — 1. 

The second line is an alternative way to specify the sum, where: 



C = {T^A\Xl5(T{dh))^l}. 



C is a subgroup of the group A of all configurations, \A\ = 2^^ in number, 
where A^ is the number of lattice sites. C is the group of all configurations 



{r} with an even number of t = —1 plaquettes occurring in the co-boundary of 
every link, i.e. forming a closed tiled surface of negative plaquettes. 
There is a second group of interest, 

C = {aeA\{a,T) = l WeC}. 

C is the group of all configurations {a} for which C{a) is different from zero. 

Therefore C is the group of a configurations which have non zero weight in 
the partition function, Eqn.(0). Further, C is the group of r configurations that 
form closed tiled surfaces as required by the explicit constraints in Eqn.fl). 

The group C has only an implicit definition here. The group C has an 
implicit definition in terms of this C. Therefore its definition is even more 
indirect. Even without an explicit definition, we have been able to specify pre- 
cisely those configurations {a} that contribute non-zero weight to the partition 
function. 

There is a third group of interest, 

V^{aeA\a = Y[ag^^}. 

where ofg^^. refers to an individual star transformation on the i'th link bi, and 
the product indicates all possible products of them. V is the group of all 
configurations {a} which can be built out of products of "star transformations" 
starting from the identity configuration. 

This is the constructive definition that is straightforward to implement in 
a simulation. We show in Ref. U that the group V is identical to the group C. 
In this way we have shown that by our proposed algorithm is ergodic. 

5 Simulation of vortex counters 

Vortices have long been considered as prime candidates for the essential dy- 
namical variable to describe confinement. A simulation offers a tool that allows 
one to correlate the occurance of vortices with values of other dynamical vari- 
ables. Hence as a first application we use this formalism to measure various 
vortex counters for Wilson loops. 

In the SO{3) x Z{2) formulation the Wilson loop is given by 

W[C] = (itr[C]775fTs)c=9S, 

where rjs and as are products of 77 and a over any spanning surface EJn'u. 

Kovacs and TomboulisEl define three vortex counters for thick vortex sheets, 
thin vortex sheets, and hybrid (patches of each on the sheet). 
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Figure 1: Fraction of Wilson loops with an even number of vortices piercing the minimal 
surface, Xe. {xo = 1 — Xe) 



Thin: 



Mhin = crs- 



If this value, = ±1, is independent of the spanning surface, then this counts 
thin vortices. 
Thick: 



M 



thick 



r]S sgn[trVK(C)]. 



This object is counting something more elusive since unlike the above case, the 
vortex structure is spread over many lattice spacings. Nevertheless it is always 
possible to find a representative of SO (3) such that the t] vortex defines the 
topological linkageH. The i] vortices can be deformed by a Z{2) transformation 
of links giving different representatives of SO (3) without cost of action. One 
can move a linked t] vortex sheet so that it no longer links the Wilson loop and 
further even transform it away. However in this case the negative contribution 
will be transferred to one of the perimeter links of the Wilson loop, and it will 
not affect the value of A/'thick- Again if this is independent of the spanning 
surface, then this counts "thick" vortices. 




Figure 2: Wilson loop, W, and tagged Wilson loop with even or odd number of "thick" 
vortices piercing the minimal surface. 



Hybrid: 



Miybrid = -^thin X -^thick = ^SVS sgn[trVF(C)]. 



As one considers all spanning surfaces, the sign of A/j^j^jj^ might change. How- 
ever if the sign of A/"^jjj|,]^ always compensates then this counts hybrid vortices. 



6 Numerical Results 

Simulations were done on a 12^ lattice for (3 ~ 2.30. Measurements were 
binned to 10 bins and jackknife errors calculated, [thin]: 200 ; [thick]: 400 ; 
[hybrid]: 200 ; [projection]: 1000 measurements. We monitor the coincidence 
of Z(2) and S0(3) monopolcs which can slip due to round off error. 

Kovacs and Tomboulis' definitions require measurements on all spanning 
surfaces. We measure here only the minimum spanning surface. Then, for 
example, A/^jiiji = — 1 does not distinguish thin from hybrid, and similarly 
■^thick " ^1 does not distinguish thick from hybrid. 

Fig. 1 shows the fraction of Wilson loops which have an even number of 
vortices, Xg, as a function of the area. All the counters approach 50% from 
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Figure 3: Wilson loop, W, and tagged Wilson loop with even or odd number of projection 
vortices piercing the minimal surface. 



above with similar behavior [xo = 1 — Xe)- They are each counting different 
things and are not expected to be equal. 

-'^hybrid i^ ^ good reference curve since it measures the sign of the Wilson 
loop itself. The area law arises from a near cancelation of fluctuating values 
due to approximately equal occurance and absence of thin, thick or hybrid 
vortices. 

The "thick" curve lies on top of the "hybrid" one. Hence the added factor 
of (Js in the hybrid counter has little effect here. We expect a vortices to be 
heavily suppressed for increasing /? since they cost action proportional to the 
vortex area. However a.t (3 = 2.30 the density of Z(2) (or S0(3)) monopoles 
= 0.2155(2) (Random plaquette signs would give a density of 0.5). Hence 
in spite of the near coincidence of these two curves, a vortices and a patches 
of hybrid vortices are important at this value of p. There is further p/idence 



below. The M 



projection 



case gives the same result as reported in RefllJ. 



Fig. 2 [Thick segment piercing the minimal surface] By definition 
We has an even number of thick segments piercing the minimal surface. Yet 
it still has an exponential fall off with area. See Fig. 4 which gives the Creutz 
ratio showing that We has about half the string tension of the W. The thin 
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Figure 4: Creutz Ratio: xi^t /) = — In — \Jtr r-Ti^ ' similarly for tagged Wilson loops 

We. 



segments are still present and they account for the disordering. 

Fig. 3 [Projection vortices piercing the minimal surface] For loops 
of area 9 and higher, Jic sign of the vortex counter correlates with Wo- These 
data agree with Rcf. i3. Comparing this with Fig. 2 [thick] the latter data are 
about a factor of 10 larger, indicating a large discrepancy in these two methods 
of tagging center vortices. 

Fig. 4 [Creutz Ratio] W and We [thick] show a constant string tension 
for larger loops. We suspect that with better statistics, We [thin] will also. We 
[hybrid] shows the vanishing of the string tension if one removes the disordering 
mechanism completely. Larger loop areas are needed to decide if the Creutz 
Ratio of We [projection] will also go to zero, or stabilize which would indicate 
that a disordering mechanism remains. 



7 Summary 

The Tomboulis formalism provides a precise way to separate the physical ef- 
fects of thin and thick vortices. Simulations in these variables allows one to 
disentangle these effects as one approaches the continuum limit where thick 
vortices are expected to dominate. 
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